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2S hyperfine structure of atomic deuterium 
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We have measured the frequency splitting between the (2S, F — 1/2) and (2S, F = 3/2) hyperfine 
sublevels in atomic deuterium by an optical differential method based on two-photon Doppler-free 
spectroscopy on a cold atomic beam. The result f^) s (2S) = 40 924 454(7) Hz is the most precise 
value for this interval to date. In comparison to the previous radio- frequency measurement we 
have improved the accuracy by the factor of three. The specific combination of hyperfine frequency 
intervals for metastable- and ground states in deuterium atom D21 = 8f^P s (2S) — /hfs(I'S') derived 
from our measurement is in a good agreement with D21 calculated from quantum-electrodynamics 
theory. 

PACS numbers: 12.20.Fv, 32.10.Fn, 32.30.Jc, 42.62.Fi 
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I. INTRODUCTION 

High-precision measurements in light atomic systems 
permit accurate tests of the quantum-electrodynamics 
theory (QED). QED calculations enter into a number of 
fundamental values related to free particles and simple 
atoms. In conventional atomic systems the accuracy of 
QED tests is restricted by insufficient knowledge of the 
nuclear structure which is the main obstacle on the way 
to improve theoretical predictions for the Lamb shift and 
the hyperfine structure in hydrogen (see e.g. 0,0)- 

The leading corrections for different energy levels com- 
ing from nuclear size effects are proportional to the 
squared value of non-relativistic wave function of the elec- 
tron at the origin: 



A-Enuci = Am ■ \ip(r = 0) 



(1) 



where the value of the coefficient An is determined 
by parameters of the nucleus and does not depend on 
atomic quantum numbers. For s-levels in hydrogen- 
like systems the non-relativistic wave function at the 
origin scales with the principal quantum number n as 
\tp(r = 0)| 2 ~ 1/n 3 . Thus, if one takes the following 
difference of the energies E(nS) corresponding to two 
different s-levels n' 3 E(n'S) — n 3 E(nS), the leading con- 
tribution of nuclear effects cancels. 

Recently, a significant progress in calculations of the 
8E KFS (2S) - E RFS {IS) difference of 25 and 15 hyperfine 
splitting (HFS) intervals in light hydrogen-like atoms has 
been achieved Q. New state-dependent QED terms to 
the HFS interval frequency /hfs (nS) up to the fractional 
order of a 4 and a 3 m e /m p have been calculated as well 
as the next-to-leading nuclear structure effects. The ac- 
curacy of this theoretical prediction for D21 exceeds now 
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the experimental accuracy which is mainly restricted by 
the uncertainty in the determination of the /hfs (25). 



Experimental results are available for the 15 and 25 
hyperfine interval in hydro gen jjj IH @, Q, deuterium 
[11} and the helium-3 ion The study of neu- 

tral atoms and ions requires different experimental tech- 
niques. While 15 HFS intervals in neutral atoms 0, Q 
have been measured with a higher relative accuracy as in 
the helium-3 ion ^(j,the 25 HFS intervals in hydrogen 
0, Q and deuterium are known less precisely as in the 
helium-3 ion [TH . Since the traditional microwave meth- 
ods have likely reached their limits we have been working 
on an optical determination of the hyperfine interval in 
the metastable 25 state of hydrogen and deuterium. 



Recently we have measured the hyperfine splitting of 
the 25 state of the hydrogen atom applying an optical 
technique • It is based on two-photon spectroscopy on 
a cold beam shielded from magnetic fields. The 25 HFS 
interval has been determined from the frequency differ- 
ence of two stable light fields exciting 15-25 transitions 
for singlet (F = 0) and triplet (F = 1) components. The 
differential method cancels a number of systematic ef- 
fects intrinsic to two-photon spectroscopy which provides 
a significant increase of accuracy in comparison with ab- 
solute frequency measurements 0, 0, [w| . 



We have improved this technique [3 and applied it 
to the spectroscopy of atomic deuterium. The hyperfine 
structure in deuterium is approximately 4 times smaller 
than in hydrogen and its accurate optical measurement is 
a harder problem, however, a smaller HFS value has also 
some advantages due to a reduction of certain system- 
atic effects for the comparison of the 15-25 frequencies 
between different spin states. The experiment and its 
systematic effects are presented in two following sections 
while a comparison of theory and experiment is summa- 
rized in the concluding part of the paper. 
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FIG. 1: Laser system for 1S-2S deuterium spectroscopy. 



II. MEASUREMENT OF THE 25" HYPERFINE 
SPLITTING 
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FIG. 2: Vacuum part of the experiment. HV and UHV are 
the high- and ultrahigh vacuum regions, while V is the bypass 
valve reducing the atomic flow escaping from the nozzle. 



The 25 HFS interval has been measured once by Reich, 
Heberle, and Kusch The value of 



fuFs( 2S ) = 40 924 439(20) Hz 



(2) 



has been obtained by a radio-frequency (rf) method ap- 
plied earlier to a similar measurement in hydrogen |5j . In 
comparison to the hydrogen measurement [5|, the abso- 
lute accuracy of determination of fj^P s (2S) 0] has been 
improved by a factor of 3. In spite of the fact that 
the HFS intervals in deuterium are strongly affected by 
magnetic fields, as in hydrogen, the average of two rf 
transitions frequencies in the metastablc deuterium atom 
(F = 1/2, m F = -1/2) (F = 3/2, in p = 1/2) and 
(F = 1/2, rn F = 1/2) <-► (F = 3/2, m F = -1/2) con- 
tains no linear field-dependent terms and thus is rather 
insensitive to the field. Still, the most important sys- 
tematic effects contributing to the 20 Hz error were the 
uncertainty in the determination of the magnetic field 
and the rf Stark effect. 

For an independent measurement of f^P s (2S) we use 
the hydrogen spectrometer setup described elsewhere 
[l5j . A dye laser operating near 486 nm is locked to 
a definite TEMoo mode of an ultra-stable cavity with 
Ultra-Low Expansion (ULE) glass spacer and a drift less 
than 0.5 Hz/s (FigQJ. To change the laser frequency with 
respect to the cavity mode we use a broadband double- 
pass acousto-optic modulator (AOM) placed between the 
laser and the cavity. The frequency of the dye laser is 
doubled in a Barium /3-borate crystal, and the result- 
ing 243 nm radiation is coupled to a linear enhancement 
cavity inside the vacuum chamber where the two-photon 
excitation takes place. 

Atomic deuterium produced in a 15 W, 2.5 GHz rf 
gas discharge flows through teflon tubes to a copper noz- 
zle cooled by a flow-through cryostat (FigJ2J|. Atoms 
thcrmalizc with the nozzle walls and leave the nozzle in 
both directions coaxially with the enhancement cavity 
mode. A part of the atoms flying along the mode are ex- 
cited to the 25* state and then reach the detection zone 
where a weak electrical quenching field is applied. The 



field mixes the 2S and 2P states of the atoms causing 
their fast decay to the ground state with the emission of 
Lyman-a photons. Photons are counted by a solar-blind 
photomultiplier connected to a photon counter incorpo- 
rated into time-of-flight measurement scheme. The exci- 
tation radiation is periodically blocked before the cav- 
ity by a chopper wheel, and 121 nm photons are de- 
tected only within a 3 ms dark time interval. Introduc- 
ing a delay r between closing the beam and starting the 
detection, we can thus select different velocity groups 
contributing to the signal. The delay r sets an upper 
limit for the atomic velocity to u max = 1/t, where I is 
the distance between the nozzle and detector. Using a 
multi-channel scaler we sort the counts in adjacent time 
bins (r = 10 fj,s, 210 jj,s, 410 ju,s, . . .,) and simultaneously 
record up to 12 spectra containing the information about 
velocity-dependent effects, e.g. 2nd order Doppler effect. 

Here we will point out some important changes which 
have been introduced in the spectrometer configuration 
(which has been used for the recent measurements in 
atomic hydrogen 0, El) to optimize it for the current 
deuterium measurement. To maximize the count rate for 
slow atoms (v ~ 200 m/s) one has to find a compromise 
between thermalization and recombination processes at 
the cold walls of the nozzle by adjusting the nozzle geom- 
etry and temperature. The best rates are observed with 
a larger nozzle diameter (2-2.5 mm compared to 1.2 mm 
for hydrogen) and at slightly higher temperatures (6-6.5 
K compared to 5 K). Typically, the count rate with the 
deuterium beam has been higher than for hydrogen (by 
a factor 2-5 depending on r). 

On the way between the gas discharge and the nozzle 
we have installed a bypass valve which can be opened to 
a high- vacuum (HV) part of the vacuum chamber, where 
the cavity mirrors, the nozzle and the detector have been 
placed. By opening the valve we reduce the atomic flow 
through the nozzle by a factor of 4. Even at the lowest 
pressure in the gas discharge (0.7 mbar) and with opened 
valve, the count rate was sufficient to detect a solid 25 
signal for the atoms with v ~ 200 m/s. The high- vacuum 
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FIG. 3: Two-photon transitions between hyperfine compo- 
nents of the IS and 25 levels in atomic deuterium. The split- 
ting of magnetic sublevels in a magnetic field is also presented. 
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volume is pumped by a turbo-molecular pump to 5 • 1CP° 
mbar. 

Atoms are collimated by two 1.5 mm diaphragms sepa- 
rating the HV and the ultrahigh- vacuum (UHV) volumes. 
Most of the interaction region (95%) is in the UHV re- 
gion pumped by a 10 4 1/s cryopump. Without the atomic 
beam, the background gas pressure in the UHV region 
equals 3 ■ 10 -8 mbar. The background gas pressure in- 
creases with the deuterium flow escaping from the nozzle 
and can be thus varied up to 2 • 10 -7 mbar. 

The excitation region is shielded against external mag- 
netic fields. As in pj, we use a two-stage /x-metal shield- 
ing. The outer shielding together with the compensation 
field reduces residual fields over the entire interaction re- 
gion to 10-20 mG. Configuration of the inner shielding 
has been improved to increase the shielding factor as well 
as its throughput for pumping. 

In the absence of magnetic fields, the 25 HFS interval 
frequency is given by the following combination of optical 
frequencies and the ground state HFS interval: 



/hfs(^) 



F=3/2 ~ /f=1/2 



(3) 



where /f=i/2 an d /f=3/2 are the frequencies of the dou- 
blet and the quadruplet transitions at 121 nm (see FigEJ. 
In our experiment, we measure the frequency difference 
(/f=3/2 — /f=i/2) using two-photon spectroscopy of the 
corresponding transitions. The frequency of the ground- 
state splitting has been measured by Wineland and Ram- 
sey || with an uncertainty of 5 • 10~ 12 : 



f^ s (lS) = 327 384 352.5222(17) Hz. 



(4) 



The contribution to the resulting error budget introduced 
by this uncertainty is negligible. 



FIG. 4: Simultaneous recording of the two transition lines 
(top). AOM frequencies /aom corresponding to the centers 
of the doublet and quadruplet transition lines (bootom) . The 
frequency change is due to the drift of the reference ULE 
cavity (bottom). 



For allowed two-photon transitions between the levels 
with AF = 0, the contribution of the linear Zeeman ef- 
fect cancels and the HFS frequency shift in a magnetic 
field B scales as -B 2 //hfs which is unfavorable for deu- 
terium in comparison to hydrogen. But due to a low-field 
regime and the absence of linear terms, our measurement 
is rather insensitive to magnetic fields. In our case, the 
magnetic field B shifts the frequency Jhfs(25) approxi- 
mately as 35 • B 2 kHz/G 2 . 

The laser as well as the lock electronic are adjusted in 
such a way, that the two (15, F = 1/2) —> (25, F = 1/2) 
(doublet) (15, F = 3/2) -> (25, F = 3/2) (quadruplet) 
transitions can be excited by changing the AOM fre- 
quency (see Fig|3J|. 

To measure the frequency difference between the two 
two-photon transitions we have applied the following pro- 
cedure. The frequency of the double-pass AOM placed 
between the reference ULE cavity and the laser is tuned 
to one of the transitions and one point of the spectrum is 
measured with a photon accumulation time of 0.5 s. After 
the measurement, the AOM frequency /aom is sweeped 
over a big fixed frequency interval (approximately equals 

to l/hfs( 15 ') - /hfs( 2S ')]/ 8 ~ 36 MHz ) to reach the 
other two-photon transition. During the sweep, the laser 

has been continuously kept in lock. The factor 8 arises 

from the double-pass scheme for the AOM, the second 

harmonic generation stage and the two-photon process. 

Than we introduce 0.5 s pause as recovery time for mul- 



4 



tiple feedbacks, and repeat the measurement of one point 
of the other two-photon transition. After the measure- 
ment we add a small frequency step for scan and sweep 
the synthesizer back to the first transition. The result of 
such a procedure is presented on FigQJtop). 

The rf synthesizer driving the AOM is continuously 
referenced to the 10 MHz signal of a HP5071A cesium 
frequency standard calibrated by GPS (FigQJ. The mis- 
calibration of the primary standard introduces a negli- 
gible uncertainty to the measured difference. After the 
double-pass AOM, the laser light is spatially filtered by 
means of a single-mode fiber, while the intensity of the 
light coupled into the cavity is actively stabilized. That 
guarantees that the coupling conditions which can change 
the frequency of the reference cavity mode |Tfjj are the 
same for both transitions. 

Both simultaneously recorded doublet and quadruplet 
transition lines are fitted with a Lorentzian function in 
the time and frequency domains. A typical time inter- 
val between the centers of such a line pair is less than 
5 s. This corresponds to some small correction on a 
level of a few hertz which should be introduced to cor- 
rect for the reference cavity drift (figQJbottom)). Such 
a procedure significantly reduces the influence of refer- 
ence cavity frequency fluctuations on the 25 HFS data 
statistics. The new method suppresses the influence com- 
ing from all drifts in the experimental setup with a time 
scale exceeding 200 s which is required to record one 
line. Even though such fluctuations do not introduce sys- 
tematic shifts, a long time is required to average them. 
The approach is much more forgiving to systematic ef- 
fects: one independent point is obtained after 3 minutes 
of measurement, while for the previous method the 
corresponding time interval was about 20 minutes. 



III. RESULTS AND SYSTEMATIC EFFECTS 

We have measured /hps (25) during 7 days using the 
new method and during 6 days applying the method de- 
scribed in 0- To test for systematic effects we have 
varied different parameters of the experiment within 
each day. We have recorded over 1000 deuterium time- 
resolved spectra. The averaged amplitude ratio between 
quadruplet and doublet transitions equals to 2.00(2) 
which is in a good agreement with theoretical expecta- 
tion. 

The systematic effects for our two-photon beam spec- 
troscopy are well-known [l3j . For the differential mea- 
surement, the most significant effects cancel out. As 
shown in J7J, the differential dynamic Stark shift can- 
cels to the level of 10~ 6 relative to the shift of the 15-25 
transition (about 500 Hz). Because of some residual fluc- 
tuations of the 243 nm radiation intensity, the doublet 
and quadruplet lines acquire slightly different shifts. We 
correct for this difference for each line by monitoring the 
power leaking out of the cavity (Fig|2J|. The resulting 
correction for the new method equals 0.5 Hz which is a 



factor 4 lower than for the measurement in hydrogen |7| . 
For the simultaneous recording technique, long-term in- 
tensity fluctuations in the enhancement cavity contribute 
less to the result than in the case of • Besides the cor- 
rection, we add a 0.5 Hz uncertainty to the error budget. 

A DC electrical field E mixes 25 and 2P levels caus- 
ing a 25 level shift. According to ^7|> the energy shift of 
the 25 level is reversely proportional to the Lamb shift. 
Thus, the differential shift scales as i?HFs/^2S-2P, where 
L2S-2P is the Lamb shift of the 25 level. The evaluated 
shift of the 25 HFS interval in deuterium equals —300 
E 2 Hz cm 2 /V 2 and is about 4 times smaller as in hydro- 
gen. The excitation region is shielded from stray fields 
by coating all surrounding parts of the interaction region 
with graphite. Residual stray fields in the setup are esti- 
mated to be less than 30 mV/cm [l8|, which corresponds 
to a 25 HFS frequency shift on a level of —0.3 Hz. We 
add an uncertainty of 0.5 Hz to the error budget con- 
sidering the slightly worse geometrical properties of our 
shielding compared to the Faraday cage used in [l8j . 

Residual magnetic fields split the magnetic sublevels 
and shift the measured hyperfine splitting to higher 
values. The shift mostly originates from imperfectly 
shielded parts of the excitation region lying in the HV re- 
gion (Fig[3J). To check for this effect we have once made a 
set of measurements without compensation field and ob- 
served a change of /^pg(25) equal to 20(22) Hz. Without 
compensation field, the measured value of magnetic fields 
around the nozzle and the detector is about 200 mG. The 
compensation field reduces the magnetic fields for at least 
a factor of 10, which means that the residual shift should 
be on a sub-hertz level. We conservatively estimate the 
shift as 0.5(1.0) Hz. 

Another source for level shifts is a pressure shift caused 
by collisions with the background gas and within the 
beam itself. The ground-state hyperfine splitting in hy- 
drogen and deuterium is rather insensitive to collisions. 
Typically, this shift is on a level of 1 Hz/mbar depending 
on the buffer gas |19l |2(j which corresponds to a vanish- 
ing shift for our pressure range. On the other side, there 
are no reliable experimental data for the pressure shift of 
the 25 HFS frequency HQ]. The upper limit for the 25 
HFS frequency shift can be taken as the frequency shift 
of the 15 (F = l,m F = ±1) ^ 25 (F' =_Lm' F = ±1) 
transition equal to -8(2) MHz /mbar ^Eg. But con- 
sidering the theoretical works j2j |^| , there is no reason 
to expect the 25 HFS shift to be orders of magnitude 
higher than for the ground state. 

We have experimentally investigated the influence of 
collisions on the fj^psi^S) frequency. Most of the colli- 
sions take place within the atomic beam where the pres- 
sure can be much higher than the background gas pres- 
sure. Using the method of simultaneous detection, we 
have varied the atomic flow over a wide range by chang- 
ing the pressure in the gas discharge (between 0.8 mbar 
and 8 mbar) and by opening the bypass valve (Fig[2J). 
The background gas pressure scales linearly with the flow 
with a small offset of 3 ■ 10~ 8 mbar. Experimental data 
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FIG. 5: The 2S HFS frequency vs. deuterium flow. One unit 
corresponds to 3.6 ■ 10 17 particles coming to the nozzle per 
second. The data is obtained using the simultaneous record- 
ing method and corresponds to the delay time of r = 410 fis. 
To extrapolate to the zero flow we fit the data with a linear 
function (solid line). 



are presented in FigJSJ Each point represents the result 
of statistical averaging of multiple data points detected 
at different days. About 1/4 of all the data have been 
taken at the lowest pressure in the gas discharge with 
an opened valve (the left point). At high flows the slow 
atoms are accelerated and pushed away from the beam 
by collisions (Zacharias effect) 18] which causes a loss 
in statistics and corresponding increase of uncertainty. 
Moreover, the measurement time becomes restricted by 
a fast growth of a film of molecular deuterium on the 
nozzle. For the previous method the flow range of {2- 
8} (Fig|5J) is practically unreachable. Disregarding the 
delay time r we have not observed any indication for a 
25 HFS interval pressure shift on our level of accuracy. 

The observed transition lines are shifted by the sec- 
ond order Doppler effect in a range between 0.1-1 kHz 
depending on the velocitydistribution of the atoms con- 
tributing to the signal [13j. The lines are not symmet- 
ric and significantly differ from Lorentzian profiles at 
short (r < 410/xs) delays (Fig^J). If we excite two- 
photon transitions in the same thermal beam and de- 
tect the signal with some precisely defined delay time 



TABLE I: 2S HFS frequency and statistical uncertainty of 
extrapolation for different delays r. The velocity « ma x repre- 
sents the maximal velocity of atoms contributing to the signal. 
The Coefficient k is the slope of the linear fit (Fig0|. 
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FIG. 6: Quadruplet transition lines for the different delays r. 



r equal for both singlet and quadruplet transitions, the 
shift and the influence of the asymmetry should cancel 
out. To estimate the possible uncertainty coming form 
the "wrong" fitting procedure with a Lorentzian function 
and all other velocity-dependent effects, we have evalu- 
ated the /hps (25) frequency for different delays. The 
data for each delay has been corrected for the AC Stark 
shift, and then the extrapolation to zero flow has been 
made (Fig[5J|. The results of this extrapolation are pre- 
sented in the Table [I] All velocity-dependent systematic 
effects, if they exist, should reveal themselves in such 
evaluation procedure. 

All the results are consistent within the error bars. For 
higher delays r the corresponding uncertainty grows due 
to the lack of statistics. As a final value we have chosen 
the result f^ ] s (2S) = 40 924453(5) Hz (r = 410 (is) due 
to the following reasons: (i) the statistical uncertainty 
for this result is small, (ii) the line-profile is practically 
undistinguishable from Lorentzian, and (iii) the fast non- 
thermalized atoms do not contribute to the signal. We 
add an additional uncertainty of 0.3 Hz corresponding to 
the conservatively estimated shift coming from the colli- 
sions with the residual background gas. 

The data obtained during 6 days of measurement ap- 
plying the older method presented in give the result of 

f$ s (2S) = 40 924462(15) Hz. The data are measured 
only at low deuterium flows {0.2-1.5} and are simply av- 
eraged over the entire ensemble. We have increased the 
uncertainty due to the possible pressure shift (1 Hz) and 
the AC Stark shift (2 Hz). 

Combining both these statistically independent results 
and adding all known systematic uncertainties linearly, 
we arrive at the final value for the 2S hyperfine interval 
in atomic deuterium 

f^> s {2S) = 40 924454(7) Hz. (5) 

This result is in a good agreement with the one obtained 
by an rf method 9], but is more accurate by a factor 3. 
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TABLE II: Summary of systematic errors and the final result 
for the 2S hyperfine interval. The independent results ob- 
tained by the "new" simultaneous recording method and the 
one from Q are presented. 
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The results as well as the error budget are collected in 
Table CU 

IV. D 2 i DIFFERENCE 

To derive the D21 difference we should combine our 
result with the known value of the 15 HFS interval 

£>£ P = 8/^(25) - f$ s (lS) = 11 280(56) Hz. (6) 

This result is in a good agreement with the theoretical 
prediction Q 

£>^ cor = 11312.5(5) Hz. (7) 



The experimental and theoretical values are presented 
in Fig[7| Restricted by the experimental accuracy, we 
can test the state-dependent QED contributions on the 
level of 0.2 ppm which can be compared to similar tests 
in hydrogen Q and the 3 He + ion We point out 

that the absolute accuracy of our measurement exceeds 
the both cited results, while the relative accuracy suffers 
from the relatively small HFS interval in deuterium. 



Today precision optical methods challenge even in 
fields were classical rf methods have been traditionally 
considered as favorable. After a successful measurement 
of the Lamb shift jH H| and the 25 the HFS interval 
in the hydrogen atom we have now presented a result 
on the 25 HFS interval in deuterium with an accuracy 
exceeding that of rf methods. 



The accuracy of the optical measurement can be im- 
proved by the further reducing the atomic velocity. The 
current results indicate that the pressure shift, which lim- 
its the accuracy of absolute frequency measurements in 
ultra-cold hydrogen ,22| , may not play such a crucial role 
for the differential technique. 
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